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ABSTRACT 

We investigate the star-formation rate (SFR) and stellar mass, M* relation of a star-forming 
(SF) galaxy sample in the XMM-LSS field to z ~ 3.0 using the near-infrared data from the 
VISTA Deep Extragalactic Observations (VIDEO) survey. Combining VIDEO with broad¬ 
band photometry, we use the SED fitting algorithm CIGALE to derive SFRs and M* and have 
adapted it to account for the full photometric redshift PDF uncertainty. Applying a SF selec¬ 
tion using the D4000 index, we find evidence for strong evolution in the normalisation of the 
SFR-M* relation out to 2 ~ 3 and a roughly constant slope of (SFR ex M“) a = 0.69 ± 0.02 
to 2 ^ 1.7. We find this increases close to unity toward 2 ^ 2.65. Alternatively, if we apply a 
colour selection, we find a distinct turnover in the SFR-M* relation between 0.7 < z < 2.0 
at the high mass end, and suggest that this is due to an increased contamination from passive 
galaxies. We find evolution of the specific SFR oc (1 + z) 2m at log 10 (M*/M 0 ) ~ 10.5, out 
to 2 < 2.4 with an observed flattening beyond z ~ 2 with increased stellar mass. Comparing 
to a range of simulations we find the analytical scaling relation approaches, that invoke 
an equilibrium model, a good fit to our data, suggesting that a continual smooth accretion 
regulated by continual outflows may be a key driver in the overall growth of SFGs. 


Key words: galaxies: evolution-galaxies: stellar content-infrared: galaxies : methods: statis¬ 
tical: infrared: galaxies 


1 INTRODUCTION 

Quantifying the interplay between star formation and stellar mass 
(M*) is a crucial component to understanding the build up of galax¬ 
ies over cosmic time. Whilst early works established the decline 
in SFR as a function of redshift (e.g. Madau 1995; Lilly et al. 
1996; Hopkins & Beacom 2006), the recent acquisition of data 
from broad band photometry coupled with access to deep and wide 
galaxy surveys, led to Noeske et al. (2007a) demonstrating a strong 
correlation between star-formation rate (SFR) and M* out to z ~ 1, 
which they termed the star forming ‘Main Sequence’ (SF-MS). 
This SF-MS was modelled by a simple power-law of the form 
SFR =f3M* where strong evolution in the normalisation results in 
a rapid decline in SFR from z ~ 1 to z ~ 0.3 (at a fixed stellar 
mass), with a slope a ~ 0.7 at z ~ 0.5. Complementary studies 
at this time by Elbaz et al. (2007) and Daddi et al. (2007) found a 
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similar trend with slopes close to unity at z ~ 1 and 2 respectively. 
Subsequent work showed an anti-correlation between the specific 
SFR and stellar mass out to z ~ 2 which has been termed down¬ 
sizing (Cowie et al. 1996) where the more massive galaxies form at 
higher redshift (see e.g. Noeske et al. 2007b; Damen et al. 2009a; 
Rodighiero et al. 2010). 

Since then, there have been many investigations of this rela¬ 
tionship, using both observations and simulations, with the aim of 
shedding light on how it connects with the underlying physical pro¬ 
cesses governing galaxy evolution. Whilst the general trends in the 
SF-MS from observations seem to agree, the way in which we ob¬ 
tain SFR measurements, coupled with the choice of initial mass 
function (IMF), the particular star-forming galaxy (SFG) selection 
criterion, the modelling of star formation histories (SFH) and stel¬ 
lar population synthesis models (SPS), dust attenuation, extinction, 
photometric redshifts, metallicities, adopted cosmology and incom¬ 
pleteness, can all introduce bias and calibration issues for mod¬ 
elling the SFR-M* relation accurately (e.g. Salimbeni et al. 2009; 
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Figure 1. Left: The VIDEO redshift distribution from our SERVS-matched 
sample of 52,812 galaxies. Right: The contribution to the redshift distri¬ 
bution by sources detected by MIPS 24 (4881 galaxies), IRAC 5.8 (909 
galaxies) and SPIRE 250 /mi (901 galaxies). 


Conroy et al. 2009; Magdis et al. 2010; Walcher et al. 2011; Pforr 
et al. 2012; Conroy 2013; Buat et al. 2014; Rodighiero et al. 2014; 
Speagle et al. 2014). 

Typical SFR indicators include the UV-wavelength range (e.g. 
Erb et al. 2006; Magdis et al. 2010; Rodighiero et al. 2011; Lee 
et al. 2011, 2012) since the emission is dominated by young mas¬ 
sive short-lived stars. The UV, however, is heavily affected by ex¬ 
tinction and so combining this indicator with IR measurements 
helps to correct for dust attenuation in the UV which is then re¬ 
emitted in the IR and FIR (e.g. Santini et al. 2009; Salmi et al. 2012; 
Whitaker et al. 2012; Reddy et al. 2012; Smith et al. 2012), or using 
the IR alone as in e.g. Oliver et al. (2010); Elbaz et al. (201 1). Neb¬ 
ula emission lines, e.g. Her, 0[ll] and 0[m] have also been used 
(e.g. Chen et al. 2009; Shim et al. 2011; Zahid et al. 2012; Suzuki 
et al. 2015), as well as from radio continuum emission and stacking 
techniques (e.g. Dunne et al. 2009; Pannella et al. 2009; Karim et al. 
2011; Zwart et al. 2014). Recently, SED modelling algorithms us¬ 
ing so-called energy balance methods that exploit the conservation 
of energy between the stellar light absorbed by dust that is then re¬ 
emitted in the mid-IR/far-IR/submm (Burgarella et al. 2005; Noll 
et al. 2009; da Cunha et al. 2008), have also been used to obtain es¬ 
timates of the total SFR as well as stellar masses and other physical 
properties. 

All of these approaches have advantages and disadvantages. 
However, a key issue is the initial selection of SFGs (e.g. Ilbert 
et al. 2010; Karim et al. 201 1). Recent studies by Rodighiero et al. 
(2014) and Speagle et al. (2014) demonstrate how different colour 
selection of SFGs leads to large variations in determining the slope 
of the SF-MS throughout the literature, making direct comparisons 
challenging. 

Nevertheless, some of the most recent work presents evidence 
for a more complex evolving MS beyond a standard power-law. 
For example, Heinis et al. (2014); Magnelli et al. (2014); Schreiber 
et al. (2015) and Whitaker et al. (2014) present results that indicate 
a MS that flattens off toward the high end of the mass function, 
favouring a MS modelled by e.g. a broken power-law. In particu¬ 
lar, Whitaker et al. (2014), using data from the 3D-HST photomet¬ 
ric catalogues and a stacking analysis of MIPS 24/xm probed the 
low mass SF-MS down to log 10 (M*/M©)= 8.4 and 9.2, at z=0.5 
and 2.5 respectively. From this analysis they parameterised the data 
with a two-power law model, and found evidence for strong redshift 
evolution for galaxies with masses log 10 (M*/M©)> 10.2. How¬ 
ever, they found no evolution for galaxies below this mass limit, 
suggesting that high-mass SFGs may follow a different evolution¬ 
ary path than lower mass galaxies. This trend at the high-mass end 



Figure 2. Comparison of the outputs from CIGALE using the full wave¬ 
length coverage available and a restricted coverage. The left-hand panels 
show a comparison of SFR estimates from CIGALE where SFRi is the re¬ 
sulting SFR using only the input photometry from VIDEO+CFHT+SERVS. 
The label SFR 2 is for the same set of objects when we include additional 
filter information from MIPS+IRA C+Herschel. The right-hand panel shows 
the same comparison for the stellar masses M* and M|. The black points 
in each panel are the medians with the ctmad errors. The grey solid 
lines on the top left and right panels indicates the one-to-one relation. The 
colour scale indicates increasing redshift from blue to red which covers 
0-1 < 2 ; < 3.0. The bottom two panels show the log 10 of the ratio between 
SFRi and SFR 2 (left) and the stellar masses (right). 


has also been reported by Tasca et al. (2014) using data from the 
VIMOS Ultra-Deep Survey (VUDS) at * ~ 1.5. 

In parallel, over the last decade there have been significant 
advances in simulations providing insight into galaxy formation 
through the modelling of physical and environmental processes. 
The most popular approaches that have emerged are the semi- 
analytical models (SAMs) and hydrodynamical simulations. SAMs 
use an analytical approach, requiring up to ~50 parameters to char¬ 
acterise the gas dynamics within and around halos (e.g. Bower et al. 
2006; Kitzbichler & White 2007; Bower et al. 2008; Guo & White 
2008). Whilst they have been largely successful at reproducing ob¬ 
served SFRs of local populations (z < 0.4), the work by e.g. Elbaz 
et al. (2007); Daddi et al. (2007); Santini et al. (2009) and Damen 
et al. (2009b) have shown they consistently under predict the star 
formation rate by a factor of 2-5 out to z ~ 2. 

With the advances in computational power, hydrodynamical 
simulations have grown in popularity. This approach attempts to 
model the physical dynamics of cosmic gas in a self consistent way 
by including the modelling of star formation processes, feedback 
from stellar winds, supernovae, and black holes (e.g. Dave et al. 
2011, 2013; Dubois et al. 2014; Vogelsberger et al. 2014; Sparre 
et al. 2014). However, while there may be broad agreement with 
overall trends of the SFR-M* relation with observations, they too 
exhibit a consistent under prediction in the overall normalisation 
across most redshifts (e.g. Daddi et al. 2007; Tasca et al. 2014). 

In this work we explore the SF-MS, focusing on using CIGALE 
(Noll et al. 2009), which incorporates empirical relations to build 
a grid of libraries that model: SFHs using SPS models (e.g. Fioc 
& Rocca-Volmerange 1997; Bruzual & Chariot 2003; Maraston 
2005), interstellar gas and dust (e.g. Calzetti et al. 1994; Calzetti 
1997; Dale & Helou 2002), fraction of AGN (e.g. Fritz et al. 2006) 
and optical and UV interstellar emission and absorption lines (see 
e.g. Kinney et al. 1996; Noll et al. 2004). We use observations pri- 
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Figure 3. Example illustrating how we propagate photometric redshift PDF 
estimates through CIGALE. The top panel shows a typical redshift PDF out¬ 
putted by lephare and the blue open circles represent the gridded redshift 
steps. The best-fit redshift for this object is shown by the vertical dashed 
line. To estimate the resulting PDF for the SFR and M*, we run CIGALE 
for a series of redshift steps that encapsulate the full range of the z p h~ PDF 
as shown in the middle and bottom panels by the red solid dots. Finally, we 
weigh the resulting SFR and M* distributions by the z p h~ PDF probabilities 
and shown in the right-hand middle and bottom panels. For completeness, 
the horizontal dashed lines denote the resulting CIGALE fitted values for the 
best-fit lephare input redshift value (i.e. a single redshift) with the corre¬ 
sponding CIGALE error estimates on these parameters shown as the vertical 
red error bars. 

marily derived from the VISTA Deep Extragalactic Observations 
(VIDEO) (Jarvis et al. 2013). 

The format of this paper is as follows. In § 2 we provide full 
details of the data and our sample selection followed by a descrip¬ 
tion of the SED fitting code CIGALE in § 3. In this section we also 
discuss how we can extend this algorithm to account for the full 
propagation of photometric redshift PDFs and discuss complete¬ 
ness estimates and other sources of uncertainties. In § 4 we discuss 
the various ways in which we model the SF-MS and then present 
our results in § 5. In § 6 we compare our observational results to a 
variety of simulations to investigate how the key physical processes 
invoked in the simulations to reproduce local relations, effect the 
SF-MS at 0.5 < z < 3. A final discussion and conclusions can 
be found in § 7. Throughout this work we adopt an Ho=70 km s -1 
Mpc -1 , Qm=03, Qa = 0.7 cosmology. 

2 SAMPLE SELECTION 
2.1 The VIDEO Survey 

For our analysis we use K s -selected galaxies from the VISTA Deep 
Extragalactic Observations (VIDEO; Jarvis et al. 2013) Survey. 


12.0 


11.5 
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Figure 4. Stellar-mass completeness as a function of redshift.The black dots 
show the distribution our final sample of 52,812 galaxies. The transpar¬ 
ent yellow region shows the stellar-mass limit M* im , for all objects that 
could be detected with an apparent magnitude limit of K s < 23.0. The 
black dotted line indicates the mass completeness limit for the K s band. 
As our sample is jointly selected with IRAC1 and IRAC2 bands at 3.6 and 
4.5 gm we compute the mass completeness limits for these bands as shown 
by the dashed and solid black lines respectively. The blue dots trace our fi¬ 
nal mass completeness limit which selects the highest mass limit from the 
three bands. 


VIDEO is a ~ 12 degree 2 survey in the near-infrared Z, Y, J, 
H and K s bands that has been designed to trace the evolution of 
galaxies and clusters as a function of both epoch and environment 
out to z rsj 4. VIDEO is being carried out over three of the most 
widely observed high-Galactic-latitude fields, covering ~3 deg 2 in 
the ELAIS-S1 field, ~4.5 deg 2 in the XMM-Newton large-scale 
structure field, and another ^4.5 deg 2 in the extended Chandra 
Deep Field-South. These fields have a wealth of data from the 
X-rays through to the radio waveband, and are, along with COS- 
MOS/UltraVISTA, the primary fields for observations with future 
facilities in the Southern hemisphere. 

In this work we focus on an area within the XMM-LSS 
field where VIDEO observations overlap with the Canada-France- 
Hawaii-Telescope Legacy Survey Deep 1 (CFHTLS-D1 Ilbert et al. 
2006) field, providing optical photometry in the ugriz bands. The 
creation of the CFHTLS/VIDEO catalogue extracted over the 10 
bands (ugriz from CFHTLS and ZJHK S from VIDEO) is de¬ 
scribed in detail by Jarvis et al. (2013). The current full catalogue 
contains 431, 949 objects over a ~ 1 deg 2 area. 

To determine photometric-redshift (z p h) estimates with 
VIDEO we used the public code LEPHARE 1 (Ilbert et al. 2006). 
For more details regarding the LEPHARE settings see Jarvis et al. 
(2013). As is relevant to this work, LEPHARE outputs a redshift 
probability distribution function (PDF) for each object, which we 
propagate in our analysis in deriving SFR and M* estimates, 

For our analysis we adopt a conservative approach and only 
select objects brighter than a limiting apparent K s -band magnitude 
K s < 23.0, which corresponds to > 95 per cent completeness for 


1 http://www.cfht.hawaii.edu/~arnouts/LEPHARE/ 
lephare.html 
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Figure 5. The two methods of selecting SFGs. Left panel, the upper plot shows the rest-frame u-r colour distribution. The bottom panel shows how this varies 
as a function of log 1 O (M*/M 0 ). The vertical dashed line indicates our colour cut to isolate SFGs; all galaxies with u-r >1.9 are considered to be passive. 
Right panel the D4000 index as determined by CIGALE. The extended tail in the D4000 distribution shows the residual passive galaxy population that we 
remove using a cut at D4000=1.3. The bottom panel shows the passive population are mostly massive objects with log lo (M*/M 0 )> 10. 
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Figure 6. Using the D4000 index to select star forming galaxies as dis¬ 
cussed in § 3.3. The colour gradient shows how the D4000 index varies as 
function of the M*-SFR distribution for a sample at z ~ 1. All parameters 
were derived from CIGALE. 


the VIDEO data (see Jarvis et al. 2013) for 0.1 < z p h < 3.0, thus 
reducing our parent sample to 74,063 galaxies. 


2.2 Matching to Multi-wavelength data 

The design of the VIDEO survey footprint was driven by the 
need to cover the wealth of ancillary data available in fields of 
several square degrees routinely accessible from Paranal. Multi¬ 
wavelength observations within VIDEO fields include Spitzer (e.g. 
Lonsdale et al. 2003; Mauduit et al. 2012) and Herschel (Oliver 
et al. 2012), in addition to ongoing ground-based optical surveys 


from the CFHT, the VLT Survey Telescope (VST 2 ) and the Dark 
Energy Survey (Flaugher 2005; Bernstein et al. 2012). For our 
analysis we thus built a well-sampled complete multi-wavelength 
dataset of galaxies detected by VIDEO and matched to available 
UV-to-IR data. 


2.2.7 SERVS Data 

The Spitzer Extragalactic Representative Volume Survey (SERVS, 
Mauduit et al. 2012) covers the three VIDEO fields and consists 
of deep observations in Spitzer- IRAC1 and IRAC2 bands at 3.6 
and 4.5 pm. We searched for SERVS counterparts to the VIDEO 
sources within a radius of 1” and only selected objects detected in 
both bands with fluxes higher than 1.0 ply (3 -a limit). Our final 
sample contains 52,812 galaxies. This jointly selected sample at 
K s , 3.6 and 4.5 pm ensures that we can reliably detect the emission 
from the old stellar populations that contribute the bulk of the stellar 
mass, up to our redshift limit. 


2.2.2 SWIRE and HerMES Data 

The Spitzer Wide-Area Extragalactic (SWIRE Lonsdale et al. 2003) 
survey provides IRAC 5.8 and 8.0 pm and MIPS 24, 70 and 160 
pm data over the whole of the VIDEO footprint. For SWIRE, we 
used the Spitzer Data Fusion (Vaccari et al. 2010 3 ) catalogues and 
matched SERVS and SWIRE positions within a search radius of 
1 ”. 

Far-infrared data was obtained from the HERschel Multi¬ 
tiered Extragalactic Survey (HerMES; Oliver et al. 2012) carried 
out with the SPIRE (Griffin et al. 2010) instrument onboard the 
Herschel Space Observatory (Pilbratt et al. 2010) at 250, 350 and 


2 http://www.mattiavaccari.net/voice/ 

3 http://www.mattiavaccari.net/df/ 
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500 /im. The SPIRE source extraction is based on forced pho¬ 
tometry at the IRAC positions of 24 ptm sources as described in 
Roseboom et al. (2010), and using an input catalogue based on the 
Spitzer Data Fusion. As with the SWIRE data we match the Rose- 
boom et al. (2010) catalogue to our VIDEO-SERYS objects within 
1 arcsec. 


3 ESTIMATING SFR AND MASS WITH CIGALE 


To estimate SFRs and M* we use the public SED fitting code, Code 
investigating GALaxy Emission, (CIGALE 4 ; Burgarella et al. 2005; 
Noll et al. 2009; Giovannoli et al. 2011), which is optimised to 
provide physical information of galaxies by fitting their UV-to-IR 
SEDs. In this section we summarise the main steps of the algorithm. 

CIGALE is based on the use of a UV-optical stellar SED plus an 
IR-emitting, dust component. In essence, CIGALE builds UV-to-IR 
models which are then fitted to the observed SEDs. The code fits the 
observed data in the UV, optical, and NIR with models generated 
with a stellar populations synthesis code, assuming a star formation 
history (SFH) and a dust attenuation as input. The energetic balance 
between dust-enshrouded stellar emission and re-emission in the IR 
is conserved by combining the UV/optical and IR SEDs. 

In our implementation, we have used the stellar population 
synthesis models of Maraston (2005) (hereafter M05) which con¬ 
sider the thermally pulsating asymptotic giant branch (TP-AGB) 
stars, and use the Kroupa (2001) initial mass function (IMF). We 
generate star formation histories (SFH) based on two stellar pop¬ 
ulations: a recent stellar population with a constant SFR on top 
of an older stellar population created with an exponentially de¬ 
clining SFR. Model spectra are then reddened using attenuation 
curves from Calzetti et al. (2000), and we use the semi-empirical 
one-parameter models of Dale & Helou (2002) (hereafter DH02) 
to fit the stellar-heated dust emission. The DH02 library is com¬ 
posed of 64 templates parametrised by a , the power-law slope of 
the dust mass over heating intensity, a is directly related to the 
iWm/fioo/im flux ratio, where feo^m and fimpm represent fluxes 
at 60 and 100 /im, respectively. CIGALE allows us to estimate the 
fraction of Lm due to an AGN by combining AGN with the star- 
burst templates. We used the set of six templates from the Fritz et al. 
(2006) library. The final part of this process uses the Kinney et al. 
(1996) and Noll et al. (2004) spectra to create empirical templates 
to consider the UV and optical interstellar emission and absorption 
lines respectively. The SFRs are computed as 


SFR = 


ATgai 

T ( e t/r _ ’ 


( 1 ) 


where M ga i is the galaxy mass, r is the e-folding (or decay) time 
and t is the look-back time. Finally, the physical parameter values 
corresponding to the best-fit model are constrained via a % 2 min¬ 
imisation. 

We explore the impact of the CIGALE outputs when flux in¬ 
formation is unavailable in certain filters for a given galaxy. The 
fraction of objects in our sample that were detected in IRAC 3 & 4, 
MIPS and SPIRE is quite small. Thus, the dominating contributing 
filters used in CIGALE are from VIDEO, CFHT and SERVS. For 
this test we have isolated those galaxies that have detections across 
the entire filter set. We then run CIGALE twice; first with input pho¬ 
tometry from only VIDEO, CFHT and SERVS, and secondly where 
we add the additional information from IRAC 3 & 4, MIPS and 


4 http://cigale.lam.fr/ 


Table 1. Input parameter for propagating the z p h PDFs through CIGALE. 
The first column shows the range of ages for the stellar population (SP) 
models used to create star formation histories. The second column shows 
the redshift range we permit based on the SP ages. 


Ages of old SP models 
r f (Gyr) 

Data z-range 
(best-fit z ph ) 

z-step coverage 
(cigale model range) 

5,6,7, 8,9 

0.1 < 2? < 0.3 

0.01 < z step < 0.41 

2.5, 3.5, 4.5, 5.5, 6.5 

0.3 < s < 0.6 

0.15 <C Zstep 's 0.75 

1,2, 3, 4, 5 

0 .6<z< 1.0 

0.45 Zstep 'n 1.15 

1,2, 3,4 

1.0 < s < 1.4 

0.85 "C Zstep "n 1.55 

1, 1.5,2, 2.5,3 

1A< z < 2.0 

1.25 Zstep 's 2.15 

1, 1.5, 2, 2.5 

2.0 < z < 2.4 

1.85 "C Zstep A- 2.55 

1.0, 1.2, 1.5, 1.7,1.9 

2A<z< 3.0 

2.25 < Zstep < 3.15 


SPIRE. Figure 2 shows the comparison between the estimated SFR 
and M* from CIGALE under these two scenarios. For the compar¬ 
ison of SFRs, we find very good agreement between the outputs, 
with an average overall scatter in the residuals of gm = 0.22 dex. 
Similarly for M* we find a very tight correlation with an average 
scatter of just gm — 0.07 dex in the residuals. We therefore con¬ 
clude that CIGALE can robustly estimate the SFR and M* without 
the need for the far-IR data to the redshift limit of our sample. 


3.1 Propagating the z p h uncertainties through CIGALE 

CIGALE requires the redshift of the galaxy as input and is not cur¬ 
rently optimised to simultaneously constrain redshift and physical 
properties of the galaxies. We therefore modified CIGALE to obtain 
both SFR and M* estimates when we have the full photometric 
redshift PDF. For our sample such information was obtained from 
LEPHARE and propagating this through the CIGALE framework re¬ 
quired only a slight modification. 

This modification is illustrated in Figure 3. The top panel 
shows a typical redshift PDF outputted from LEPHARE, where the 
blue circles represent the resolution of the LEPHARE grid for a 
given galaxy. As we can see in this example, there are two peaks 
in the z p h~ PDF; the first at 0.3 and the second at z~0.6 with 
a lower probability. The majority of the contribution is centred at 
z~0.3 and LEPHARE returns a best-fit value as indicated by the 
vertical dashed line. To include the full z p h contribution for a given 
galaxy we run CIGALE at a series of redshift steps that cover the 
range that fully describes the z p h- PDF. The middle and bottom 
left-hand panels in Figure 3 show the resulting log 10 (SFR) and 
log 10 (M*) estimates from CIGALE respectively as a function of 
redshift. Finally, we can weight this relation in SFR and M* ac¬ 
cording to the interpolated probability distribution of the z p h , as 
shown in the right-hand bottom and middle panels. We also show 
the estimated error on the SFR and M* from CIGALE for this galaxy 
using the best-fit redshift value. 

To account for the range in the full photometric redshift PDF 
estimation from LEPHARE, we ran CIGALE with different stellar 
population (SP) model ranges according to a set of redshift slices 
which are summarised in Table 1. The use of applying two stel¬ 
lar populations was explored in more detail recently by Buat et al. 
(2014), who applied a similar approach to derive SFRs and M* es¬ 
timates at redshifts z > 1. Our approach allows for z p ^-PDFs that 
may not be well described by e.g a simple gaussian. In § 5 we apply 
this approach in determining SF-MS. 
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redshift 

Figure 7. Cosmic variance as a function of redshift for a range of mass bins 
for the current VIDEO data release covering 1 deg 2 . 

3.2 Stellar mass completeness 

To estimate the completeness limit for M* we follow the approach 
used by Ilbert et al. (2013) and Pozzetti et al. (2010), and which 
is illustrated in Figure 4, where we show the stellar masses for our 
sample of 52,812 galaxies as a function of redshift. We calculate 
the lowest stellar mass which could be detected for a given galaxy 
with 

log 10 (M, im ) = log 10 (M.) + 0.4 (K s - (2) 

where, in our study, the apparent magnitude limit. A - ' 1 " 1 = 23.0. 
This is shown as the transparent yellow region on Figure 4. We then 
bin the data in redshift and compute the 90th percentile of stellar 
mass within each bin, shown as blue dots in the figure. Thus for 
a given redshift, we can interpolate the stellar mass completeness 
limit that corresponds to the mass for which 90 per cent of galaxies 
fall below the percentile. 

3.3 Selecting star-forming galaxies 

To analyse the SF-MS we need to isolate SFGs and remove any 
passive galaxies that would contaminate the sample. However, se¬ 
lecting SFGs from multi-wavelength data can be done in a variety 
of ways. This was recently reviewed by Speagle et al. (2014) and 
explored in detail by Rodighiero et al. (2014), where they high¬ 
lighted how different colour selections can lead to different SFG 
populations being isolated, contributing to the broad range of de¬ 
rived slopes in the literature. For example, applying a BzK selec¬ 
tion (e.g. Daddi et al. 2004; Rodighiero et al. 2011; Kashino et al. 
2013) tends to isolate more active SFGs, biasing against galaxies 
which would otherwise be classified as SF by other methods (e.g. 
NUVr J). This type of selection generally leads to steeper derived 
slopes for the SF-MS. To help avoid selecting only the bluest most 
active galaxies, Ilbert et al. (2013) used a mixed approach and ap¬ 
plied two colour selections in r — J and NUV-r, which provided a 
broader range of SFGs. In a recent paper by Renzini & Peng (2015) 
they suggest a construction of a 3D SFR-Mass-Number distribution 
may help in defining a more objective SF-MS by adopting a cut at 
the ridge line of the star-forming peak. 

In this study we explored two approaches which yielded very 


different results in our final analysis. With the filter coverage avail¬ 
able to us we found a rest-frame u-r was the optimal colour selec¬ 
tion approach, but not necessarily the most effective way to remove 
the passive population. The top-left panel of Figure 5 shows the 
u-r distribution where we observe a bimodal distribution, split ap¬ 
proximately at u-r ~ 1.9, as indicated by the vertical dashed line. 
The bottom-left panel shows the u-r as function of stellar mass (as 
determined using CIGALE) which perhaps better illustrates this seg¬ 
regation in the two populations. As such we removed all objects in 
our final sample with u-r > 1.9. However, as we shall discuss in 
the § 5, this does lead to a significant fraction of passive galaxies 
being missed, along with potential SFGs being removed from the 
sample at redshifts z <1. 

Alternatively, one may use the 4000 A break (hereafter re¬ 
ferred to as D4000). The D4000 index provides useful informa¬ 
tion regarding the age of the stellar population and perhaps pro¬ 
vides a more physical selection process. A low index indicates a 
young stellar population, whereas the larger the value the older the 
galaxies (e.g. Kauffmann et al. 2003; Westra et al. 2010). We use 
CIGALE to determine the D4000 index as part of the fitting pro¬ 
cedure, adopting the definition as detailed in Balogh et al. (1999). 
Moreover, the calculation is intrinsically a dust-free measurement, 
as it is based on the unreddened spectra. Whilst this approach is 
SED-model dependant, it does provide us with a self-consistent ap¬ 
proach to identify SFGs. In the right-hand panel of Figure 5 we 
show the resulting D4000 distribution as a function of M*. We ap¬ 
ply a conservative cut at D4000 =1.3, and consider galaxies beyond 
this limit to be older passive galaxies. For the remainder of this arti¬ 
cle our results will be based on this D4000 selection. As we can see 
in Figure 6 this type of selection imposes a much sharper cut in the 
SFR-M* plane as a function of redshift impacting on conclusions 
drawn about the evolution of the SF-MS as we will discuss in more 
detail in § 5.4. 


4 MODELLING THE SF-MS 

We parameterise and model the SF-MS in a number of ways to test 
the robustness of our findings and to allow direct comparison with a 
number of related works. The most common approach is to assume 
a simple power-law of the form 

lo Sio( SFR ) = a lo Sio( M *) + Pi ( 3 ) 

(see e.g. Noeske et al. 2007a; Daddi et al. 2007; Elbaz et al. 
2007; Santini et al. 2009; Heinis et al. 2014). Without a redshift 
dependent term, the data can be simply split into a series of redshift 
slices, and a fit is perfomed within each slice to determine how the 
normalisation [3 and the slope a of the SF-MS evolve. 

We calculate the total error on the median SFRs in each bin 
as cr to t —\/(^mad + a v)’ where ctmad is the error on the median 
(median absolute deviation, or MAD), and is related to the stan¬ 
dard error on the mean a by a < x > = ax/1.4826 (Hoaglin et al. 
1983). This is combined in quadrature with the cosmic variance 
uncertainty cry. This is the uncertainty in observed number density 
of galaxies arising from the underlying large-scale density fluctua¬ 
tions. We have used the public code GETCV developed by Moster 
et al. (2011), to estimate cosmic variance as a function of redshift 
for the range in M* considered. Figure 7 shows how the fractional 
cosmic variance error cry varies as a function of redshift for a range 
of stellar mass bins as computed for our sample. We then perform 
a least-squares fit to estimate the power-law parameters in Equa¬ 
tion 3. 
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Figure 8. The SF-MS for SFGs from 0.1 < z < 3.0. The coloured regions are a 2D histogram of the SFR-M* distribution, where the number of objects 
in each pixel (dSFR = 5M*= 0.2) is indicated in the colour bar at the top of the figure The white filled circles denote the median values with their 

respective uncertainties computed as at ot=y ^(^mad + which combines the mean absolute deviation and cosmic variance (see text for details). The 
grey area shows the region of our stellar mass completeness limit. Any data within this region is not considered in the final analysis. iVg a i indicates the total 
number of galaxies in each panel, with the total number within our stellar mass completeness range shown in parenthesis. The correlation coefficient between 
SFR and M* for the stellar mass-complete data is shown as r c in each panel. 
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Figure 9. The dispersion in the SF-MS, jms, as a function of redshift for 
a series of stellar-mass bins. The grey shaded region shows a range of ctms 
as reported in e.g. Noeske et al. (2007b); Whitaker et al. (2012) and Mag- 
nelli et al. (2014). The region between the dashed line indicates our average 
dispersions across all stellar mass bins where for log 10 (M*/ Mq)=10.0, 
10.5, 11.0 and 11.5, (Jms= 0.29, 0.29, 0.31 and 0.28 respectively as esti¬ 
mated within the stellar mass completeness limits. 

Although using the median SFR can be useful when dealing 
with a statistically-small sample, it naturally leads to a loss of in¬ 
formation. Therefore, we also perform power-law fits within each 
redshift slice to the full galaxy distribution which we refer to as 
‘all-data’ fits and consider the uncertainty on the SFR and M* for 
each individual galaxy as determined using CIGALE. 

Recent studies by e.g. Magnelli et al. (2014); Whitaker et al. 
(2014); Schreiber et al. (2015) and Tasca et al. (2014) have found 
a flattening of the MS towards the high mass end, and therefore 
require more robust modelling using either e.g. a broken power-law 
or an extra term in the form of a quadratic given by 

logio(SFR) = Ai + A 2 log 10 (M„) + A 3 [log 10 (M*)] 2 . (4) 

Others have instead chosen to model the SF-MS jointly with a 
redshift-evolution term. In Schreiber et al. (2015) they parameterise 
as 

lo Sio[ SFR ( z )] = lo g 10 (M») - ai + /3i log 10 (l + z) 

- /3 2 [max(0, log 10 (M») - a 3 log 10 (l + z)] 2 (5) 

Alternatively, Whitaker et al. (2012), applied a second order poly¬ 
nomial of the form 

logi 0 [SFR(z)] = a(z)[log 10 (M*) - 10.5] + /3(z), (6) 

where the terms a(z), /3(z) are given by, 
a(z ) = ai + 0 L 2 Z, and 
P{z) — Pi + ($ 2 % + 03 z 2 , 

and a±, a. 2 , Pi, P 2 and P 3 are the fitting parameters. To explore 
evolution we adopt the Whitaker et al. (2012) form, applying our 
z p h- PDF approach as detailed in § 3.1. In practice we create 1000 
Monte Carlo (MC) samples randomly drawn from these PDFs, 
which are then individually fitted via Monte Carlo Markov Chain 
(MCMC). We can then examine the distribution of parameter con¬ 
straints to assess any bias from the standard fitting procedure. We 
have used the public MCMC algorithm CosmoMC (Lewis & Bri¬ 
dle 2002). Samples are generated using the Metropolis-Hastings 


algorithm (Hastings 1970) for each of several intercommunicating 
Markov chains, with a ‘proposal distribution function’ that allows 
full exploration of the posterior space. We adopt a x 2 likelihood 
given by 

- 2togl0 (£, sx ‘-± N:.ii( SFR .‘““> ~ loSnISFRi '--)] 1 


where log 10 (SFR data ) is the observed data and log 10 (SFR™ od ) 
is the model. We apply this model to our z p h ~PDF analysis of 
the SF-MS and our cr to t in this case includes the intrinsic disper¬ 
sion of the SFR-z and M*-z distributions as well as cr v giving 

<7tot (z) = 0) 2 + CTsfr(^) 2 + CTy]. 

4.1 Calibration of the main-sequence 

In § 5.2 we compare our SF-MS results to a range of recent studies. 
However, as we can see from Table 2 there is a spread in the under¬ 
lying assumptions to determining stellar masses and SFRs that in¬ 
clude: wavelength coverage/selection of an observed sample, SFR 
estimation (L — SFR relation), the IMF, stellar population syn¬ 
thesis models (SPS), extinction, metallicities, adopted cosmology, 
not to mention other contributing factors such as the modelling star 
formation histories (SFH), dust attenuation, photometric redshifts 
and incompleteness. Whilst for any given study the results may 
be self consistent, quantifying differences in the SF-MS beyond 
their intrinsic scatter due to differing assumptions can be challeng¬ 
ing. In general, calibration of MS data is typically an average cor¬ 
rection which impacts the normalisation. In an attempt to make a 
fair comparison we have adopted a similar calibration approach de¬ 
tailed in Speagle et al. (2014) (SI4). In their study they examined 
twenty-five recent studies and highlight three major sources that 
require a calibration offset, namely the choice of IMF, SFR indi¬ 
cators (L — SFR relation) and SPS modelling. For the works that 
overlap with ours in Table 2 we have used the S14 re-calibrated 
MS power-law fits which we show in Table 4. Whilst S14 estimate 
offsets for e.g. adopted cosmology, extinction, emission line, these 
were shown to be sub-dominant to the IMF, SPS models and SFR 
indicators. As such, in our calibration we consider only contribu¬ 
tions from the latter. 

The IMF offsets used in S14 were taken from taken from Za- 
hid et al. (2012) where they calibrated to the Kroupa IMF using 

Mf = 1.06M? = 0.62 M?, (8) 

where the superscripts refer to Kroupa, Chabrier (Chabrier 2003), 
and Salpeter (Salpeter 1955) IMFs, respectively. As they calibrate 
to the Kroupa IMF, these would lead to respective offsets in stellar 
mass for Chabrier and Salpeter of = +0.03 and Cm* = -0.21. 
This would translate to an offset in the normalisation of the MS fit 
by 

Pc — Po — ol x Cm, , (9) 

where Pc is the calibrated normalisation, Po is the original fitted 
normalisation and a is the slope. Since we have applied the Kroupa 
IMF in our analysis we do not include this offset in our calibration. 

There have been a number of investigations into the differ¬ 
ing approaches to generating SPS models and in particular the 
treatment of the modelling of thermally pulsating asymptotic gi¬ 
ant branch (TP-AGB) stars, (e.g. Salimbeni et al. 2009; Conroy 
et al. 2009; Magdis et al. 2010; Walcher et al. 2011). To calibrate 
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Figure 10. Parameter constraints from our MCMC analysis of the model shown in Equation 6. left panel set: The dark and light blue shaded contours show 
the 2-D lcr and 2o marginalised results respectively. The red and yellow contour regions show an example from the 1000 Monte Carlo z p h ~PDF constraints. 
The 1-D marginalised probability distributions are shown on the diagonal as the blue line (red line) for the VIDEO median ( 2 ^-PDF) fits. The cross near 
the centre of each contour set shows the best-fit parameter value, right panel set: Comparison of Whitaker et al. (2012) median fits (red and yellow shaded 
contours) with our median fits (blue contours). 


from M05 to BC03 (as in S14) we adopt the conservative offset of 
Ci£ s = +0.15 dex which was an estimated average from analyses 
from Magdis et al. (2010); Salimbeni et al. (2009); Conroy et al. 
(2009). Therefore, we apply Equation 9 to obtain the MS offset. 

In a recent study by Buat et al. (2014) they examined poten¬ 
tial SFR offsets derived from CIGALE compared to that of a stan¬ 
dard L — SFR relation that convert FUV and IR luminosities in to 
SFRs, and found a very tight correlation with that of the CIGALE 
derived SFRs. In the case of the two population model, which we 
have applied in our analysis, they reported a small offset resulting in 
a slight underestimate of the instantaneous SFR with SFRir 5 fuv- 
For our analysis we therefore do not attempt a correction, but intend 
to investigate this more deeply in subsequent work. 

To summarise, we have found that our dominant source of cal¬ 
ibration required for our comparison is from the SPS models. Our 
total calibration offsets in the normalisation of our ‘raw’ MS fits 
can be found on Table 4.1 (shown in parenthesis) along with the 
S14 derived offsets for the works we compare to in Table 2. 


5 RESULTS 

Figure 8 shows the SF-MS for our SFG sample across twelve red- 
shift bins. The grey region in each panel shows the M* complete¬ 
ness limit; galaxies within this region are excluded from our anal¬ 
ysis. The blue and white circles show the median SFRs in bins of 
width AM* = 0.2 dex. 

We find a consistent correlation between SFR and M* out to 
2 : < 3.0 which we quantify with the correlation coefficient, r c , 
in each redshift panel on the Figure. We find an average value 
of r c — 0.61; a trend consistent with many other works (e.g. 

Erb et al. 2006; Daddi et al. 2007; Noeske et al. 2007b; Dunne 
et al. 2009; Pannella et al. 2009; Rodighiero et al. 2010; Elbaz 


et al. 2011; Whitaker et al. 2012; Muzzin et al. 2013; Magnelli 
et al. 2014). We have also estimated the dispersion of the SF- 
MS denoted as cjms in Figure 9, for four stellar mass bins with 
Alog 10 (M*/M o ) = 0.5 dex, centred on log 10 (M*/M©)= 10.0 
(red), 10.5 (blue), 11.0 (green) and 11.5 (black). Previous studies 
by e.g. Noeske et al. (2007b); Whitaker et al. (2012) and Magnelli 
et al. (2014) have reported an average dispersion over all stellar 
mass ranges between ~0.3 and 0.35, which is indicated by the 
shaded grey region in Figure 9. As we can see, our values for cfms 
vary about this region as a function of redshift, where in general at 
z < 1.4 we find 0.27 < a ms < 0.32, at 1.4 < 2 < 1.6 the scatter 
increases to cfms ~ 0.38 before showing an overall decline towards 
z ~ 3 for the highest stellar mass bins with log 10 (M*/M©)=11.0 
and 11.5. This decline in dispersion was also alluded to in Whitaker 
et al. (2012) and Daddi et al. (2007). If we take the average values 
of cfms for each stellar mass bin across their respective redshift 
ranges we find for log 10 (M*/M©)=10.0, 10.5, 11.0 and 11.5 dis¬ 
persions of aMS= 0.29, 0.29, 0.31 and 0.28 respectively. 


5.1 Main Sequence constraints -1 

In this section we model the SF-MS relation using the redshift- 
dependent model given by Equation 6 and to check the robustness 
of our results we model our data with three approaches. In the first 
we explore the contribution from our z p h -PDF derived values for 
SFR and M* as detailed in § 3.1. We then model our observed 
full sample that includes intrinsic errors from CIGALE and finally 
fit to the median values of Figure 8 which allows a more direct 
comparison to the work of W12. Table 3 summarises the results of 
these fits. 

In the second column in Table 3 we account for each LEPHARE 
z p h -PDF by creating a set of 1,000 Monte Carlo (MC) data-sets by 
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Table 2. Details of the properties of the various constraints we compared to in Figure 12. Column 2: the initial mass function (IMF) used, column 3 is 
the indicator used estimate SFRs, column 4: wavelength selection of the parent sample, column 5: method to select star-forming galaxies, column 6: stellar 


population synthesis (SPS) model used, column 7: extinction law, column 8: adopted cosmology, column 9: survey and column 10: total area of sample. 

Source IMF SFR Indicator Wavelength SFG SP Model Extinction Cosmology Survey Area 

selection selection Law (h, ) deg 2 

This work 

Kroupa 

SED (CIGALE) 

K/ 3.6, 4.5 /Ltm 

D4000, u-r 

M05 

d C00 

(0.7,0.3,0.7) 

VIDEO/CFHT/SERV S 

1 

Daddi et al. (2007) 

Salpeter 

Ha 

K 

sBzK 

BC03 

e P84 

(0.7,0.3,0.7) 

GOODS 

0.033 

Dunne et al. (2009) 

Salpeter 

1.4 GHz 

K 

None 

BC03 

COO 

(0.71,0.27,0.73) 

UDS/UKIDSS 

0.8 

Elbaz et al. (2007) 

Salpeter 

Ha 

UV/opt 

blue 

a PEGASE.2 

f CF00 

(0.7,0.3,0.7) 

SDSS/GOODS 

705 

Elbaz etal. (2011) 

Salpeter 

FUV 

FIR 

blue 

PEGASE.2 

COO 

(0.7,0.3,0.7) 

GOALS/AKARI/SDSS 

38960 

Heinis et al. (2014) 

Chabrier 

UV+IR 

UV 

None 

BC03 

None 

(0.7,0.3,0.7) 

HerMES/COSMOS 

1.68 

Magnelli et al. (2014) 

Chabrier 

UV+IR 

K 

3<r-clip 

BC03 

None 

(0.71,0.27,0.73) 

COSMOS 

1.97 

Noeske et al. (2007b) 

Kroupa 

UV+IR 

Optical 

redshift 

BC03 

9 CE01 

(0.7,0.3,0.7) 

AEGIS 

1 

Oliver et al. (2010) 

Salpeter 

IR 

FIR 

None 

b R08 

None 

(0.7,0.3,0.7) 

SWIRE 

11.33 

Pannella et al. (2009) 

Salpeter 

1.4 GHz 

sBzK 

sBzK 

BC03 

P84 

(0.7,0.3,0.7) 

COSMOS 

0.9 

Reddy et al. (2012) 

Salpeter 

UV+IR 

UV 

LBG 

c CBll 

COO 

(0.7,0.3,0.7) 

various 

0.081 

Rodighiero et al. (2010) 

Salpeter 

UV+IR 

4.5 /urn 

U-B 

BC03 

COO 

(0.7,0.27,0.73) 

GOODS 

? 

Rodighiero et al. (201 1) 

Salpeter 

UV+IR 

NIR+FIR 

sBzK 

BC03 

COO 

(0.73,0.26,0.74) 

COSMOS/GOODS-MUSIC 

1.73 

Santini et al. (2009) 

Salpeter 

UV+IR 

zlKI 4.5 ^m 

2cr-clip 

M05 

C97/C00 

(0.7,0.3,0.7) 

GOODS-MUSIC 

0.04 

Schreiber et al. (2015) 

Salpeter 

UV+IR 

H/K 

UVJ 

BC03 

COO 

(0.7,0.3,0.7) 

CANDELS-GOODS/COSMOS 

1.8 

Whitaker et al. (2012) 

Kroupa 

UV+IR 

NIR 

UVJ 

BC03 

COO 

(0.7,0.3,0.7) 

NMBS 

0.4 

Whitaker et al. (2014) 

Chabrier 

UV+IR/H/3 

NIR 

UVJ 

BC03 

COO 

(0.7,0.3,0.7) 

CANDELS/3D-HST 

0.25 

Zahid et al. (2012) 

Chabrier 

Ha/H/3 

Optical/UV 

Lines 

BC03 

COO 

(0.7,0.3,0.7) 

SDSS-DR7/DEEP2 

0.5 


a Fioc & Rocca-Volmerange ( 1 997), b Rowan-Robinson et al. (2008), c Charlot & Bruzual (2011) (never formally published) 
d Calzetti et al. (2000), e Prevot et al. ( 1984), f Chariot & Fall (2000), 5 Chary & Elbaz (2001) 


Table 3. SFR-MS parameterised according to Equation 6. The second column shows the results when we propagate the full redshift PDFs through CIGALE, 
computed as the mean of the 1000 Monte Carlos fits with their respective averaged ltr standard deviations. The third column is fitted to the full observed 
data-set including the intrinsic errors on stellar mass and SFR from CIGALE (cr? ig ), the uncertainty on the best-fit redshift from lephare (<rf eph ), the intrinsic 
dispersion of the SFR and stellar mass distributions (c^ isp ), and finally the uncertainty due to cosmic variance cry 2 . The forth column shows our fits to the 
median values and the last column shows the fits to the median values from Whitaker et al. (2012). The reduced x 2 /V is shown in the bottom row. 

This Work Whitaker et al. (2012) 

(0.1 < z ph < 3.0) (0 < z < 2.5) 

Parameter '-*-' 


(Equ. 6) 

z p h ~PDF (all data) 

[ CT ? ig + ± ph +< T disp +< T V 2 ] 

(medians) 

(medians) 

a l 

0.68 db 0.01 

0.68 ± 0.01 

0.62 ± 0.04 

0.73 ±0.16 

OL2 

0.004 ±0.011 

0.004 ±0.011 

0.03 ± 0.04 

-0.12 ±0.16 

01 

0.33 ± 0.02 

0.29 ± 0.02 

0.24 ± 0.05 

0.41 ±0.17 

02 

1.36 ±0.03 

1.45 ±0.03 

1.47 ±0.08 

1.11 ±0.28 

03 

-0.30 ± 0.01 

-0.32 ± 0.01 

-0.33 ± 0.03 

-0.18 ±0.11 

x 2 A 

0.70 

0.78 

0.49 

0.05 


randomly sampling from the PDFs and generating corresponding 
SFRs and stellar masses as described in § 3.1. We then fit Equa¬ 
tion 6 to each MC sample where the likelihood term includes an 
additional uncertainty for the intrinsic dispersions of SFR-z and 
M*-z distributions as well as the contribution from cosmic vari¬ 
ance cr v . In this way we can directly inspect the distribution of 
constrained parameters and estimate their confidence limits. Here, 
the marginalised parameter constraints are the mean of the 1000 
MC sample distributions and the associated errors are the mean la 
standard deviations. In the third column of the table we fit Equa¬ 
tion 6 to our actual observed sample to incorporate statistical un¬ 
certainties from the intrinsic errors from CIGALE on SFR and M* 
the uncertainty on the best-fit redshift from LEPHARE, the intrin¬ 
sic dispersion of the SFR-z and M*-z distributions, and the uncer¬ 
tainty due to cosmic variance. Our median results are then shown 
in the fourth column with the W12 constraints presented in the last 
column of the table. It should be noted that we have reanalysed 
their median data points using the MCMC routine to provide error 
estimates on the parameters that were not previously available. 

In the left hand panel of Figure 10 we show an example of the 


2-D la and 2a marginalised results of one of the MC z p h ~PDF fits 
compared to our full data set and median constraints. We find that 
they are consistent to within 2a. On closer inspection of the z p h~ 
PDFs outputted from LEPHARE we found only a small fraction of 
objects in our final sample were poorly constrained and therefore 
their contribution to this fit would have minimal impact. Neverthe¬ 
less, this remains an important quality check of our sample to the 
robustness of the final result. Furthermore, comparing our median 
results to those of W12 (right panel set of Figure 10) and Figure 11, 
we find our results are consistent to within their la errors, although 
we constrain the parameters much more strongly due to the larger 
sample. 

5.2 Main Sequence constraints - II 

Table 5 shows our resulting power-law constraints as applied to the 
redshift bins in Figure 8. The term ‘all-data’ refers to our fits within 
each redshift slice to the full mass-complete galaxy distribution and 
the values in parenthesis are our re-calibrated MS fits. Firstly we 
note that the median and ‘all-data’ fits are fully consistent within 
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Figure 11. SFR versus redshift for the SFGs in our sample, in stellar mass 
bins of \og 10 (M*/Mo) = 9.5, 10.0, 10.5 and 11.0. In each panel the green 
points show the galaxy distribution, and the blue circles denote the median 
values used in the fit from Figure 8. The redshift-dependent fit to the median 
points (Equation 6) is shown as the blue solid line. The black dashed line is 
the best fit from Whitaker et al. (2012) shown for comparison. 

la and, as also found in the previous section, we observe strong 
evolution in the normalisation (/3) consistent with the picture of 
an overall increase in star formation rate, at a given stellar mass, 
with redshift. In particular we observe a steady decline in SFR for 
a fixed stellar mass of log 10 (M*/M©)~10.0 by a factor of ~ 19 
from z ~ 2 to 2 ; ~ 0.2. In the following sections we compare 
our results to the broader literature that coincide at three selected 
redshift ranges 0.1 < z < 0.8, 0.8 < z < 1.2 and 1.9 < z < 2.1. 
We use the calibrated MS fits as discussed in § 4.1 which are shown 
in Table 4 and Figure 12. In all three panels our results are shown 
as the blue line with the error on the fits shown as a shaded blue 
region. 

5.2.7 Comparison at z ~ 0.45 

In Figure 12 we show our fits to the ‘all-data’ case (denoted by 
the blue shaded region). In the left-hand panel of the figure our 
data covers a redshift range of 0.10 < z < 0.80 with a stellar 
mass complete sample of 7487 galaxies. At this redshift we find 
[a, /3] = [0.70 =t 0.01,-6.46 =t 0.06], which are fully consistent 
with Noeske et al. (2007a) who found [0.67 =L 0.08, —6.00 =L 0.78] 
(black solid line), Santini et al. (2009) [0.70 ± 0.14, -6.27 ± 1.54] 
(dashed black line) and Oliver et al. (2010) [0.73±0.05,—6.86± 
0.55] (red squares). The radio stacking results by Dunne et al. 
(2009) find a much steeper gradient to ours with constraints 
[a,/3] = [0.83 =b 0.01,-7.83 zb 0.13]. In the case of Whitaker 
et al. (2012), their data was fitted using the redshift-dependant poly¬ 
nomial model of Equation 6, which is represented by the yellow 
solid line. Since their redshift binning is not aligned with ours we 
cannot make a straightforward power-law fit comparison, however, 
as we can see in the left panel of Figure 12, their fit seems very 
closely aligned with that of Noeske et al. (2007a) in both normal¬ 
isation and slope over this redshift range. In Figure 13 we show 
the power-law fit to their data for each redshift bin (shown as black 
solid circles) and we find good agreement with our result over this 
redshift range. 


Table 4. Power-law contraints from Figure 12. The numbers in parenthesis 
are those calibrated by Speagle et al. (2014) and which were used in the 
plotting of the Figure. Details of our calibration can be found in § 4.1. 


Redshift 

a 

0 

0.1 < z < 0.8 

This work 

0.70 ± 0.01 

-6.36 (-6.46) ± 0.06 

Dunne et al. (2009) 

0.83 db 0.01 

-7.92 (-7.83) ±0.13 

Noeske et al. (2007b) 

0.67 =b 0.08 

-6.19 (-6.00) ±0.78 

Oliver et al. (2010) 

0.73 ± 0.05 

-7.02 (-6.86) ± 0.55 

Santini et al. (2009) 

0.70 ±0.14 

-6.33 (-6.27) ± 1.54 

0.8 < 2 ; < 1.2 

This work 

0.70 =b 0.01 

-5.92 (-6.02) ± 0.07 

Dunne et al. (2009) 

0.83 ± 0.05 

-7.38 (-7.39) ± 0.53 

Elbaz et al. (2007) 

0.90 

-8.14 (-8.06) 

1.9 < z < 2.1 

This work 

0.84 ± 0.04 

-7.01 (-7.13) ±0.28 

Daddi et al. (2007) 

0.90 

-7.60 (-7.52) 

Dunne et al. (2009) 

0.88 =b 0.03 

-7.21(-7.36) ± 0.37 

Pannella et al. (2009) 

0.95 db 0.07 

-8.30 (-8.35) 

Reddy et al. (2012) 

0.97 zb 0.05 

-8.28 (-8.32) ± 1.28 

Rodighiero et al. (2011) 

0.79 

-6.42 (-6.36) 

Santini et al. (2009) 

0.85 ±0.17 

-7.24 (-7.15) ± 1.87 

Zahid et al. (2012) 

0.46 db 0.07 

-2.99 (-2.94) ± 0.70 


In a similar way, Magnelli et al. (2014) modelled their data 
with a quadratic given by Equation 4. In this case we have averaged 
their results over their two redshift ranges 0.20 < z < 0.50 and 
0.50 < z < 0.80 to obtain an approximate comparison. Whilst 
their data supports a turnover toward the higher stellar masses, 
there appears to be a good agreement up to log 10 (M*/M©)~10.5 
at which point their fit diverges. Although we note that no er¬ 
rors on their parameter constraints were reported in their analysis. 
Schreiber et al. (2015) also model their data in such a way as to 
allow for a SFR turnover (Equation 5) which we can observe at this 
redshift similar to that of Magnelli et al. (2014). 

As with Whitaker et al. (2012), the binning of Rodighiero et al. 
(2010) are also not aligned with ours. However, we show our fits to 
their data in Figure 13, and for illustrative purposes have interpo¬ 
lated their fits to show their results in Figure 12. As a cross-check 
we fitted their data to the redshift-dependent model Equation 6 and 
found these fits to be consistent with our interpolated results. At this 
redshift their results show a slightly shallower slope of a = 0.60 
and normalisation a = —5.38. 

It is perhaps not too surprising that there is, in general, very 
good agreement between the results shown here at this relatively 
low redshift, since there is less likely to be residual contamination 
from passive galaxies which could result in shallower slopes. 

5.2.2 Comparison at z ~ 1 

As we move to the intermediate redshift bin 0.80 < z < 1.2 
shown in the middle panel of Figure 12, we have a total of 7746 
galaxies in our stellar mass complete sample. Here we find very lit¬ 
tle difference from our result at z ~ 0.45 where [a, /3\ = [0.70 zb 
0.01, —6.02zb0.07]. However, we now begin to see more divergent 
trend between data-sets. For example, Elbaz et al. (2007) (black 
solid line) found a slope a = 0.9, close to unity and a significantly 
lower normalisation of f3 = —8.14. A similar trend was found by 
Dunne et al. (2009) with [a,/3] = [0.83 zb 0.05, -7.39 zb 0.53] 
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Figure 12. Comparison of the SF-MS. Our fits to the data are shown by the thin blue line with the error shown as the light blue shaded region. For each panel 
we find a slope a = 0.70 =b 0.01, 0.70 zb 0.01 and 0.83 ± 0.02 at z ~ 0.45, 1 and 2 respectively. Where possible we have shown comparative measurements 
within their respective stellar mass ranges from Noeske et al. (2007b); Daddi et al. (2007); Elbaz et al. (2007); Santini et al. (2009); Dunne et al. (2009); 
Pannella et al. (2009); Rodighiero et al. (2010); Oliver et al. (2010); Elbaz et al. (2011); Rodighiero et al. (2011); Zahid et al. (2012); Reddy et al. (2012); 
Whitaker et al. (2012); Heinis et al. (2014); Magnelli et al. (2014) and Schreiber et al. (2015). It should be noted that the symbols are not data points but merely 
indicate the power-law fit. The fits by Heinis et al. (2014) in the middle panel were estimated at 2 ; ~ 1.4 using a broken power-law (see text for details). 


Table 5. Parameter constraints for log 10 (SFR) = a log 10 (M*) ± f$. The redshift slices correspond to those in Figure 8. Columns 2 and 3 the fits to the 
median values of SFR. Columns 4 and 5 are fitted to the full distribution of data points. The values in parenthesis are our re-calibrated fits from M05 to BC03 
SPS models (see § 4.1 for details). 


median fit all data fit 







-N 

/- 


-\ 

Redshift slice 


P 


a 

p 


a 

0 

1 — 1 

0 

A 

is* 

IA 

0.30 

- 7.86 (- 

-7.98) 

±2.30 

0.82 ± 0.24 

-7.64 (-7.76) 

± 0.20 

0.80 ± 0.02 

0.30 < * < 

0.50 

-6.27 (- 

-6.37) 

±2.18 

0.68 ± 0.22 

-6.35 (-6-45) 

±0.13 

0.69 ±0.01 

VI 

V 

0 

10 

0 

0.70 

-5.65 (- 

-5.74) 

±2.35 

0.63 ±0.23 

-5.73 (—fj-82) 

±0.13 

0.64 ±0.01 

0 

A 

0 

A 

is* 

IA 

0.90 

-4.79 (- 

-4.88) 

±2.05 

0.56 ± 0.20 

-5.60 (-5-69) 

± 0.11 

0.65 ±0.01 

0 

Jo 

0 

A 

is* 

IA 

1.10 

- 6.10 (- 

- 6 . 21 ) 

±2.23 

0.71 ±0.21 

-5.69 (- 5 - 79 ) 

± 0.11 

0.68 ± 0.01 

VI 

V 

0 

T —1 

1.30 

-5.75 (- 

-5.86) 

±2.59 

0.69 ± 0.25 

-5.86 (-5.96) 

± 0.12 

0.71 ±0.01 

VI 

V 

0 

CO 

i—l 

1.50 

-4.94 (- 

-5.03) 

±3.49 

0.63 ±0.33 

-5.07 (-5.16) 

±0.17 

0.64 ±0.02 

1.50 < * < 

1.70 

-4.94 (- 

-5.03) 

±3.79 

0.64 ± 0.35 

-5.35 (-5-45) 

±0.18 

0.68 ± 0.02 

1.70 < s < 

1.90 

-5.47 (- 

-5.57) 

±4.26 

0.69 ± 0.39 

-5.93 (-6.04) 

±0.23 

0.73 ± 0.02 

1.90 < * < 

2.10 

-5.93 (- 

-6.04) 

±4.53 

0.73 ±0.41 

cT 

1 — 1 

iV 

0 

0 

1 

±0.28 

0.83 ± 0.03 

2.10 < s < 

2.30 

-7.36 (- 

-7.49) 

±5.30 

0.87 ±0.48 

-7.86 (-8.00) 

±0.40 

0.92 ± 0.04 

2.30 < * < 

3.00 

-7.44 (- 

-7.58) 

±4.94 

0.88 ± 0.44 

-6.54 (- 6 . 66 ) 

±0.33 

0.80 ±0.03 


(green triangles). In contrast, Rodighiero et al. (2010) find a shal¬ 
lowing slope of a to 0.52. At this redshift our results are more 
aligned with Whitaker et al. (2012) and Magnelli et al. (2014). We 
also note that Santini et al. (2009) found a slope that becomes shal¬ 
lower between 0.6 < z < 1.5 from a = 0.73 at 2 ~ 0.8 to 
a = 0.65 at z ~ 1.2. With our D4000 SFG selection criteria we 
do not find strong evidence for a turnover in the SF-MS relation, 
however, we find that over this redshift range our results are con¬ 
sistent with those of Magnelli et al. (2014) (black dotted line) and 
Schreiber et al. (2015) (right triangles), at least up to M*~ 10.5. 
The results shown by Heinis et al. (2014) are consistent with ours 
with a modelled power-law giving [—5.7±0.7,0.70±0.07]. How¬ 
ever, they found a broken power-law a better fit to their data at this 
redshift finding a slope of ~0.85 for M*< 10 lo ‘ 5 M© and a turn¬ 
ing over to ~0.5 at the high mass end (shown as yellow squares 
in the Figure). In Section 5.4 we will explore the sensitivity of our 
modelling with respect to how SFGs are selected. 


5.2.3 Comparison at z ^ 2 

Finally, in the third panel we show the SF-MS over the redshift 
range 1.90 < z < 2.10, with 1683 galaxies considered in the 
fit. Here we find significant positive evolution of the slope towards 
unity with [a, /3] =[0.84 ± 0.04,—7.13 ± 0.28]. Here we find sepa¬ 
rations in the magnitude of the slopes. In general Rodighiero et al. 
(2010), Whitaker et al. (2012), Zahid et al. (2012) and Magnelli 
et al. (2014) favour a much shallower slope of a ~ 0.45, whereas 
we, along with Daddi et al. (2007); Dunne et al. (2009); Santini 
et al. (2009); Rodighiero et al. (2011); Pannella et al. (2009) and 
Reddy et al. (2012), find the opposite, with a range 0.79 < a < 
0.97. Schreiber et al. (2015) also reported a steeper relation towards 
higher redshifts with a turnover only observed in the lower redshift 
regime. It should be noted that the uncertainties on the parame¬ 
ters for the power-law fits within each redshift bin for the Whitaker 
et al. (2012) data are large, of the order a a ~ 0.75, which means 
that their results are formally consistent with ours at this redshift. 
The differences in a are perhaps better seen in Figure 13, which 
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Figure 13. Comparison of the power-law slope a constraints as a function 
of redshift. The blue line and shaded error region is our fiducial result when 
we apply our 794000 selection. The legend indicates comparative works 
from Noeske et al. (2007b); Daddi et al. (2007); Elbaz et al. (2007); Santini 
et al. (2009); Dunne et al. (2009); Pannella et al. (2009); Rodighiero et al. 
(2010); Oliver et al. (2010); Elbaz et al. (2011); Rodighiero et al. (2011); 
Zahid et al. (2012); Reddy et al. (2012), Whitaker et al. (2012) and Whitaker 
et al. (2014) as also compared in Figure 12. The Whitaker et al. (2014) data 
represent fits for log stellar masses < 10.2 (solid red diamonds) and > 10.2 
(solid black diamonds). 


Table 6. Parameter constraints for (sSFR /yr x ) = A + 7 log 10 (1 + z) for 
a range of fixed stellar masses as shown in Figure 14. 


Stellar mass range Redshift range A 7 

loglo(M*/Mo) 


9.85 < M* 

< 

10.15 

0.25 

< 

2 

< 

1.35 

-10.05 ± 0.01 

3.56 ± 0.01 

10.15 < M* 

< 

10.45 

0.27 < 

z 

< 

1.95 

-10.05 ± 0.01 

3.09 ± 0.03 

10.45 < M* 

< 

10.65 

0.36 < 

z 

< 

2.44 

-10.03 ± 0.01 

2.60 ± 0.04 

10.65 < M* 

< 

10.85 

0.65 

< 

z 

< 

2.86 

-9.88 ± 0.01 

2.15 ± 0.04 

10.85 < M* 

< 

11.15 

0.86 

< 

z 

< 

2.96 

-9.94 ± 0.01 

2.13 ± 0.06 


shows that the magnitude of a in our sample clearly begins to rise 
at z > 1.5 and e.g. Whitaker et al. (2012) begins to fall. 

5.3 Evolution of specific star-formation rate 

An alternative way to probe star formation histories over cosmic 
time is to look instead at the specific star formation rate (sSFR). 
This is simply given by sSFR=SFR/M* which allows us to probe 
star formation per unit stellar mass. In Figure 14 we show how 
the sSFR evolves as a function of redshift for our sample of SFGs 
for stellar-mass bins, centred on log 1 O (M*/M 0 )= 10.0,10.3,10.5, 
10.75 and 11.0. The redshift ranges we probe within each mass bin 
are consistent with our mass completeness limits. We again com¬ 
pare our findings with other available sSFR measurements from the 
literature, as indicated in Figure 14. Our results are entirely consis¬ 
tent with other work out to z ~ 3 within our stellar-mass complete¬ 
ness limits, with our data showing a flattening off in the evolution 
of the sSFR beyond z > 1.5 for our highest stellar-mass bins. 

We parameterise the redshift evolution of the sSFR for each 
stellar mass bin with (sSFR / yr -1 ) = A + 7 log 10 (l + z). These 
results are shown in Table 6 . We find that between mass bins 10.3 
and 10.5 the relation evolves from 7 = 3.09 to 2.60 out to a 
redshift of 1.95 and 2.44, respectively. This result is consistent 
with that of the recent work from Tasca et al. (2014) who found 
7 s= 2.8 =t 0.2 over a similar redshift range for fixed stellar mass 



redshift 


Figure 14. SSFR evolution. On each panel we show our median values of 
SSFR as a function of redshift for a range of log stellar masses centred 
on log lo (M*/M 0 )=lO.O,lO. 3 ,lO. 55 ,lO .75 and 11.0, shown as black, red, 
blue, yellow and green shaded regions respectively. The statistical uncer¬ 
tainty is the combined error on the median and the cosmic variance i.e. 

crtot=^j (^mad + ( 7 v ) 1 The other sSFR measurements are from Feulner 
et al. (2005); Noeske et al. (2007b); Daddi et al. (2007); Dunne et al. (2009); 
Kajisawa et al. (2009); Rodighiero et al. (2010); Karim et al. (2011); Salmi 
et al. (2012); Whitaker et al. (2012); Ilbert et al. (2013); Zwart et al. (2014); 
Tasca et al. (2014), and Schreiber et al. (2015). 


oflog 10 (M*/M©)~ 10.3 using YUDS. Overall, we find a general 
flattening of the relation with redshift as we move to our highest 
stellar mass bin of log 10 (M*/M©)=11 with 7 = 2.13 ± 0.06 and 
normalisation of A = —9.94 =L 0.01 between 0.85 < z < 2.96. 

If we restrict the redshift range out to z < 1.4, following Ilbert 
et al. (2014), we find a similar, but slightly steeper evolution with 
7=[3.39±0.01, 3.57±0.04, 3.39±0.05, 3.58±0.22, 3.91±0.36] for 
respective mass bins of log 10 (M*/M©)= 9.5-10.0, 10-10.5, 10.5- 
11 and 11-11.5. In their work they found a similar mass-dependent 
evolution in sSFR with 7 =[2.88±0.12, 3.31±0.10, 3.52=L0.15, 
3.78±0.60]. Why we find steeper slope at log 10 (M*/M o )= 10 re¬ 
mains unclear. However, we do note that at this mass range our fit 
was limited to z < 1.35 due to our completeness limits. If we relax 
our SFG D4000 selection cut from e.g. 1.3 to 1.35, which would 
introduce an older population of galaxies into our sample, we find 
that this has only marginal affect. Alternatively, if we instead use 
a u-r colour selection as discussed in § 3.3, this only really affects 
the high mass bins which is consistent with our findings in § 5.4. 
Of course there remains the possibility of numerous options within 
CIGALE that could be explored in more detail, but we reserve this 
for future work. 

Nevertheless, within the redshift range we are considering in 
this work, we find very good agreement across a broad range of re¬ 
lated works, that shows an in increasing mass-dependent evolution 


MNRAS 000, 1-19(2011) 

































14 Johnston R, et al. 


with stellar mass below z ~ 2, which then begins to slow consid¬ 
erably beyond this redshift. 

5.4 Star-forming selection effects and the high-mass turn 
over 

As we discussed in § 4.1, the majority of off-sets added to cal¬ 
ibrate the MS relation changes only the normalisation. However, 
the slope of the SF-MS is sensitive to effects such as how SFGs are 
selected (e.g. Oliver et al. 2010; Whitaker et al. 2012; Rodighiero 
et al. 2014; Speagle et al. 2014; Sobral et al. 2014; Ilbert et al. 
2014). As shown in Figure 13 there is a broad range in the con¬ 
strained values of a as a function of redshift across the literature. 
In the left panel of Figure 15 we compare our D4000 cut (red re¬ 
gion) with the rest-frame u-r colour selection (blue region) and we 
can see that between 0.4 < z < 1 there is a distinct shallowing of 
a applying our colour cut, which is a result of the flattening of the 
SF-MS at stellar masses >10.5 M©. 

This is shown more clearly in Figure 16 where we show the 
SFR versus M* relation, over the redshift slice 0.8 < z < 1.2. The 
left-hand panel shows the full SFR-M* distribution prior to any 
selection for SF galaxies. In the middle panel we show the resulting 
distribution once we impose the u-r colour cut, and finally the right 
panel shows our D4000 cut. It is clear that with a u-r colour cut 
there remains a residual population that, with the D4000 selection, 
would be considered a more passive population. This results in the 
median points having lower SFRs toward the high mass end. Thus, 
the resulting fit is shallower than that of the D4000 cut. Using the u- 
r we thus find similar model constraints to those of Magnelli et al. 
(2014) and Whitaker et al. (2014) applying the quadratic model 
of Equation 4. This fit is shown as the blue line where we find 
[A1,A2, A3] = [—24.596,4.626, —0.206]. Over the same redshift 
range Magnelli et al. (2014) found [-26.51, 4.77, -0.202] and in 
Whitaker etal. (2014) they found [-27.4±1.9, 5.0 ±0.4, -0.22 ± 
0.02] at 0.5 < * < 1.0 and [-26.0 ± 1.7,4.6 ± 0.3,-0.19 ± 
0.02] at 1.0 < z < 1.5. Although our results using a u-r colour 
selection show consistency in the model parameter fits to previous 
work, this is more likely due to residual contamination by an older 
sub-population of galaxies due to applying the u-r cut. 

In addition, if we apply a more relaxed cut in D4000 from 
1.3 to 1.35, thus allowing a slightly older galaxy population into 
our sample we can also observe the turnover in our median trends. 
This is demonstrated in the right hand panel of Figure 15, where the 
black stars are our fiducial result of D4000< 1.3 and the red squares 
are from Whitaker et al. (2014). As we can see, with the cut at 
D4000 <1.3, our results are consistent with a constant slope out to 
high mass. However, with a cut at D4000 <1.35, the slope flattens 
at log 10 (M*/M©)>10.5. These tests highlight our sensitivity of 
our fits to the initial selection of SFGs. 

Lastly, if we look at Table 2 the majority of works we compare 
to (including our own) have selected samples over an area of < 
1 deg 2 . At these scales uncertainties due to cosmic variance can 
become important particularly at the high mass end. In our work 
and that of e.g. Ilbert et al. (2013) this has been included in the 
model constraints, however, it is unclear as to whether other works 
consider this effect in their analysis. 

6 COMPARISON WITH SIMULATIONS 

In this section we compare our results for the SF-MS with a variety 
of simulations, performed with a number of different methods, to 


help elucidate the key physical processes that dictate the evolution 
of massive galaxies. 

6.1 Hydrodynamical simulations 

From the hydrodynamical approach we look at the Horizon (Dubois 
et al. 2014), Illustris (Sparre et al. 2014) and the Dave et al. 
(2013) (D13) simulations. As well as hydrodynamics and gravity, 
this genre of simulations typically include varying prescriptions of 
star formation, gas cooling and heating, and feedback from stel¬ 
lar winds, supernovae and AGN. Discussion on the relative merits 
of the various simulations is beyond the scope of this paper, how¬ 
ever, we do note that they each probe varying volumes and resolu¬ 
tions: Illustris use a box size of 106.5 Mpc /i _1 on the side with 
1820 3 DM particles and 1820 3 gas cells with particle and gas tar¬ 
get masses of 6.26 x 10 6 M© respectivley. Horizon have similar 
volume of 100 Mpc h~ 1 on the side with 1024 3 dark matter (DM) 
particles of mass 8 x 10 7 M©, and D13 used the Gadget-2 N- 
body + Smoothed Particle Hydrodynamics Springel et al. (2005) to 
evolve their simulations, within a volume of 32 Mpc h~ x on the 
side containing 512 3 gas particles of mass 4.5 x 10 6 M©. 

6.2 Scaling relation approaches 

Behroozi et al. (2013) parameterised the stellar mass-halo mass 
(M*-Mh\ SMHM) scaling relation to model observed star- 
formation histories out to z — 8. By applying a MCMC approach to 
constrain their 15 free parameters, they begin with an initial guess 
for the SMHM relation, then find galaxy growth histories by apply¬ 
ing the SMHM relation to dark matter merger trees. They then de¬ 
rive a series of inferred stellar mass functions and SFRs onto which 
they apply simulated observational errors and biases. Finally they 
compare their model to available data to constrain the SMHM rela¬ 
tion. 

The recent work by Mitra et al. (2014) (Ml4) is an analyti¬ 
cal based approach which also exploits the use of observed scaling 
relations. In their approach they begin with a mass balance equa¬ 
tion to model the inflows and outflows in the interstellar medium 
(ISM) assuming galaxies grow along a slowly-evolving equilib¬ 
rium between accretion, feedback and star formation (e.g. Finlator 
& Dave 2008). In essence, the parameters of the equilibrium model 
describe the motion of gas into and out of galaxies, which is termed 
the baryon cycle. In their modelling they have only 8 free param¬ 
eters that are constrained via a MCMC, comparing their model to 
the M*-Mh constraints of Behroozi et al. (2013) and Moster et al. 
(2013), as well as the M*-metallicity scaling relation from e.g. 
Andrews & Martini (2013); Zahid et al. (2014) and Steidel et al. 
(2014), out to z ~ 2. 

Thus, by design, these types of approaches should be able to 
better reproduce the observed SFR-M* relations compared to hy¬ 
drodynamical and SAM simulations. 

6.3 Results 

In Figure 17 we show the SF-MS at z = 1 (left) and 2 (right). The 
blue solid line shows the fit to our data above our stellar mass com¬ 
pleteness limit. The dotted blue line is an extrapolation of this fit 
for illustration purposes. The simulations show general agreement 
with the slope of the relation, but all have a lower normalisation in 
SFR by varying degrees. Perhaps as expected, the scaling relation 
approaches are most aligned with our work. By averaging over the 
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Figure 15. Left: Comparison of derived power-law slopes (a) as a function of redshift for our different SFG selection. The blue is for the D4000 j 1.3 cut 
as was shown in Figure 13 and the red shows the results for a rest-frame u-r selection. Right: Comparing the high-mass turnover. The black stars show our 
fiducial result using an SFG criterion of D4000 < 1.30. The blue circles show how our results change if we allow a slightly higher cut at 1.35, allowing older 
galaxies to be including in the SF-MS. We compare this to Whitaker et al. (2014) data shown as red squares which presented results with flattening of the 
SF-MS log 10 (M*/M o ) > 10.2 dex. 



log(M* [M 0 ]) 


Figure 16. Comparison of the impact of different SFG selection methods on the SF-MS. The left panel shows the full sample at z ~ 1 that includes a passive 
population of galaxies. In the middle panel we show a star-forming sample applying a rest-frame u-r colour cut. The white solid circles are the median points 
of the distribution. The red solid line shows the resulting power-law fit for the stellar mass-complete sample yielding a slope of a = 0.29; the dashed line 
is just the extrapolation of this fit below the completeness limit. The blue line shows a quadratic fit to the same median points which yielded coefficients of 
Al, A2, A3=[-21.527, 4.063, -0.1807] (see Equation 4). The right-hand panel now shows the star-forming sample if we use a D4000 index cut. The red line 
here shows a power-law slope of a = 0.66. 


stellar mass bins at 10, 10.5 and 11.0 we find M14 and B13 to be 
lower by a factor ~ 1.8. In contrast we find that the SF-MS from the 
hydrodynamics simulations of D13 and Illustris to be lower by as 
much as ~ 3 with Horizon showing the largest discrepancy being 
lower by a factor of ~6. At z = 2 there is improvement in the nor¬ 
malisation of M14 and B13, with their models now being lower by 
factors of 1.4 and 1.2 respectively between M* of 10.5 and 11. The 
hydrodynamics predictions show only a slight change from z = 1, 
and continue to significant under-predict the normalisation in the 
SF-MS relation. 

It is worth noting that the Horizon simulation data set available 
to us, contains both SFGs and passive galaxies. Although we are 
unable to segregate them explicitly, we artificially impose a cut in 
SFR that closely mirrors our D4000 selection, shown in the bottom 
panels of Figure 17. As we can see there is only a marginal upward 


shift in the normalisation at z — 1 and very little change at z — 2, 
implying that the bulk of the simulated galaxies have significantly 
lower SFRs at a given stellar mass compared to the observations. 

We now compare the observational results on the evolution of 
the sSFR (Figure 18), to the simulations, for the same range of stel¬ 
lar mass bins considered in § 5.3. If we consider redshifts 2 < 1.0 
across all stellar mass bins, we find that M14, B13 and Illustris 
show good agreement with our results. At 1 < z < 2 the Illustris 
predictions show an average factor of ~ 3 lower normalisation at 
log 10 (M*/M©)=10, which seems to reduce to 2 towards the high¬ 
est stellar mass bin M*= 11. We find both M14 and B13 continue to 
show very good agreement with our results across all stellar mass 
bins at z > 1. Finally, we again note that the Horizon simulations 
are lower by a factor of ~ 6 across all mass bins. 
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Figure 17. Comparison of the SFR-MS with current simulations. The shaded distribution of points with contours show our SFR-Mass results at z ~ 1 (left) 
and z ~ 2 (right). The fit to our data is shown as the solid blue line, with the dashed line extrapolated below our stellar mass completeness limit. The results 
from simulations at these redshifts are Illustris (Sparre et al. 2014, red solid line, hydro), Mitra et al. (2014) (green solid line, equilibrium model), Horizon 
(Dubois et al. 2014, black diamonds, hydro), Dave et al. (2013) (green dots, hydro) and Behroozi et al. (2013) (yellow triangles). 


7 DISCUSSION 

Theoretical predictions from the ACDM modelling of the average 
growth rate of dark matter halos through smooth cold accretion was 
well characterised to be oc (1 + z) 2 ' 25 at z rsj 2 (Dekel et al. 
2009). Developments in this area now favour a stronger evolution 
in sSFR to oc (1 + z) 2 ' 5 (Dekel et al. 2013; Faucher-Giguere et al. 
2011 ) which while lower than our constraints, is beginning to show 
convergence toward a higher evolution. In our comparison of the 
SFR-M* and sSFR-z trends with the hydrodynamical simulations 
of Illustris (Sparre et al. 2014), Horizon (Dubois et al. 2014) and 
Dave et al. (2013) we found in all cases an under prediction of the 
normalisation by a factor ranging between ~ 2 (llustris and Dave 
et al. 2013) and 6 (Horizon). Such a disparity between observations 
and simulations, both in hydrodynamical and semi analytical mod¬ 
elling (SAMs), has been identified and discussed extensively (e.g. 
Daddi et al. 2007; Elbaz et al. 2007; Santini et al. 2009; Damen 
et al. 2009b; Dave et al. 2013; Sparre et al. 2014; Genel et al. 2014; 
Tasca et al. 2014). The general consensus shows at low redshifts 
(z < 0.5) simulations can reproduce the sSFR and again at z > 4. 
However in the intermediate regime there remains a disconnect in 
the evolution of sSFR. Whilst there are a number of possible issues 
that may account for this effect that have been suggested such as: 
oversimplified gas accretion modelling, systematic offsets in gas 


cooling rates, insufficient sub-grid models that control star forma¬ 
tion and stellar feedback, this remains an unresolved issue. More¬ 
over, hydrodynamical and SAMs typically require ~ 20 and ~ 30- 
50 parameters respectively which can make drawing physical con¬ 
clusions challenging. 

Fundamentally, observations seem to be increasingly favour¬ 
ing the downsizing scenario, whereas simulated star formation his¬ 
tories are tied to the dark matter accretion histories and cannot pro¬ 
duce the observed mass-dependence in sSFR (Damen et al. 2009b; 
Somerville & Dave 2014; Sparre et al. 2014). 

The scaling relation based approaches of Behroozi et al. 
(2013) and Mitra et al. (2014) proved they could recover well the 
sSFR evolution across all mass ranges. This is perhaps not too sur¬ 
prising since their models are being constrained by current obser¬ 
vational data. Nevertheless, they do present the opportunity to ex¬ 
plore key physical aspects of galaxy formation within a relatively 
simple analytical framework compared to the full blown architec¬ 
ture of hydro and SAMs. In particular, the equilibrium model of 
Mitra et al. (2014) does not explicitly model halos, cooling, merg¬ 
ers or a disk star formation law, but instead parameterises the mo¬ 
tion of gas into and out of galaxies, assuming that the gas reservoir 
in a galaxy is slowly evolving with an equilibrium between accre¬ 
tion, feedback, and star formation (Finlator & Dave 2008; Bouche 
et al. 2010). Their results suggest that mergers are sub-dominant 
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for overall galaxy growth and that the key drivers for the average 
evolution of SFGs is the continual smooth accretion regulated by 
continual outflows with occasional cycling. 

Of course the way in which we select our SFGs plays a crucial 
role in the derived slopes of the SF-MS. In this work we demon¬ 
strated how two different SFG selection criteria can lead to very 
different constraints in the modelled SF-MS and, in particular, the 
high-mass turnover. From our comparisons to a wide range of re¬ 
lated works there is a clear and well known bi-modal distribution in 
the slope of the SF-MS particularly beyond z ~ 1; where some find 
a around 0.5 and others close to unity. In Whitaker et al. (2014) 
using UV+IR derived SFRs they found a non evolving low mass 
(log lo (M*/M 0 )< 10.2) MS slope close to unity. In contrast, they 
found an evolving and shallower high mass (log[M*/M 0 ] >10.2) 
slope that evolved to a = 0.59 at z ~ 2.25. A similar turn-over 
has also recently been reported by Heinis et al. (2014); Magnelli 
et al. (2014); Schreiber et al. (2015). In Tasca et al. (2014) this 
has been interpreted as evidence of a more gradual star formation 
quenching progressing from high masses at high redshifts toward 
low masses at low redshifts. In the recent work by Buat et al. (2014) 
showed that the type of SFR indicator affects the normalisation of 
the MS and not the slope. However as discussed in e.g. Speagle 
et al. (2014); Ilbert et al. (2013) they way in which we select SFG 
can impact the slope considerably. 

Our main approach was to select SFGs based on the D4000 in¬ 
dex as outputted by CIGALE. This selection imposed a rather sharp 
cut in the SFR-M* distributions which led to a roughly constant 
slope being derived out to z ~ 1.7 and little evidence for a high 
mass-turn over. However, when we instead applied a rest-frame u-r 
colour cut to isolate SFGs, we found strong evidence for a turn over 
at log 1O (M*/M 0 )~ 10.5 which was consistent with the trends of 
the aforementioned works. However, we note that in both scenarios 
we found slopes of a ~ 0.8 beyond z > 2 that showed no sign of 
flattening at high mass which seems more inline with the trends of 
Schreiber et al. (2015). 

The issue of which colour selection one should adopt to ob¬ 
tain a clean SF-MS has also been highlighted in e.g. Ilbert et al. 
(2013) and Speagle et al. (2014). Whilst using a mixture of differ¬ 
ent colour cuts to remove purely quiescent galaxies as in e.g. Ilbert 
et al. (2013) ( r-J, NUV-r) does seem to provide a cleaner selection, 
the use of a D4000 cut within the SED fitting framework potentially 
gives us a more physical way of identifying an ageing population. 
This has been explored in recent work by Morokuma-Matsui et al. 
(2015) where they found the D4000 break as a cleaner way to se¬ 
lect SFGs, albeit at low redshifts z < 0.2. Clearly this remains an 
ongoing issue that requires further investigation. 


8 CONCLUSIONS 

In this study we have used CFHTLS visible ugriz and VIDEO 
near-infrared Z, Y, J , H, K s photometry, along with the IRAC 
3.6 and 4.5 pin i bands from SERVS to investigate the evolution in 
the relation between stellar mass and star-formation rate to z ~ 3 
using the public SED fitting code CIGALE to determine the SFR 
and stellar mass. Where available we included observations from 
SWIRE and HerMES for the longer wavelength data. CIGALE en¬ 
forces an energetic balance between dust-enshrouded UV stellar 
emission and its re-emission in the IR, making the most of our 
multi-wavelength dataset. 

We summarise our findings as follows: 

• We propagated the full photometric redshift uncertainties in 



redshift 


Figure 18. Same as Figure 14 but now comparing to simulations. Predic¬ 
tions are from Behroozi et al. (2013), Mitra et al. (2014), (Illustris, Genel 
et al. 2014; Sparre et al. 2014) and (Horizon, Dubois et al. 2014) 


the form of their PDFs through CIGALE to provide more robust 
measurement of the SFR and stellar masses. In this work, the PDFs 
were estimated using LEPHARE and in general our results showed 
good consistency with an analysis which does not consider the 
PDFs, with no evidence for residual bias. 

• By modelling the data with a simple power law we have found 
strong evolution in the normalisation of the SFR-M* relation out to 
z ~ 3. For a fixed stellar mass of log 1 O (M*/M 0 )=10.0 we find the 
median SFR decreases by a factor of ~ 19 from z~2toz~0.2. 
This is consistent with the picture of quenching of star formation 
and downsizing. 

• By selecting SFGs using a D4000 index cut we observe a 
nearly constant slope for the SF-MS out to z ~ 1.7 of a = 
0.69±0.02, beyond which we observe arise, with a = 0.90±0.03 
at z ~ 2.2. At the highest redshift range probed, 2.3 < z < 3.0 we 
find slightly lower values for the slope of a = 0.80 d= 0.03. 

• However, when selecting on the basis of rest-frame u-r colour 
rather than the D4000 index, we found a slope that became progres¬ 
sively shallower out to z < 1.7, reaching a low value of /3 ~ 0.35 
at z ~ 1. We found that in this scenario the SFR-M* relation was 
more appropriately modelled by a broken power law or a quadratic, 
showing a turnover at log lo (M*/M 0 )~lO.5. This effect was also 
observed by relaxing our D4000 cut from 1.3 to 1.35 and thus al¬ 
lowing a slightly older population into our SF sample. On closer 
inspection we found that the u-r selection is not as efficient at re¬ 
moving the passive population and the turnover is largely due to a 
residual quiescent population that is not isolated. However, beyond 
z > 1.7 our results are consistent with the trends of the D4000< 1.3 
selection, showing a steepening of a ~ 0.8 out to z < 3. 

• Our study of the sSFR with redshift from z = 0.1 to 3 found 
strong evolution. By modelling the data as sSFR oc (1 Yz) 1 we find 
a similar mass-dependent evolution out to z < 1.4 with 7 = 3.39 zb 
0.01 to 7 = 3.91 =b 0.36 from log lo (M*/M 0 )=9.75 to 11.25 re- 
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spectively. Extending this analysis to higher redshifts yielded 7 = 
2.60 =b 0.04 for a stellar fixed mass of log 10 (M*/M©)=10.5 out to 
z ~ 2.44. At our highest stellar mass bin of log 10 (M*/M©)=ll 
we found a flattening in the relation with 7 = 2.13 zb 0.06 between 
0.86 < * < 3. 

Finally, our current sample size has been limited to 1 deg 2 to 
allow comprehensive broadband coverage. However, the VIDEO 
survey will have a total coverage of 12 deg 2 on its completion, 
allowing us to probe a larger and deeper area and thus will provide 
some of the most concise measurements of the SFR-M* relation, 
leading to greater insight into how this relates to the underlying 
dark matter distribution and its evolution with redshift. 
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